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ABSTRACT 

One of the most exciting near-term prospects in physics is the potential discovery of 
gravitational waves by the advanced LIGO and Virgo detectors. To maximise both 
the confidence of the detection and the science return, it is essential to identify an 
electromagnetic counterpart. This is not trivial, as the events are expected to be 
poorly localised, particularly in the near-term, with error regions covering hundreds 
or even thousands of square degrees. In this paper we discuss the prospects for finding 
an X-ray counterpart to a gravitational wave trigger with the Swift X-ray Telescope, 
using the assumption that the trigger is caused by a binary neutron star merger 
which also produces a short gamma-ray burst. We show that it is beneficial to target 
galaxies within the GW error region, highlighting the need for substantially complete 
galaxy catalogues out to distances of 300 Mpc. We also show that nearby, on-axis 
short GRBs are either extremely rare, or are systematically less luminous than those 
detected to date. We consider the prospects for detecting afterglow emission from an 
an off-axis GRB which triggered the GW facilities, finding that the detectability, and 
the best time to look, are strongly dependent on the characteristics of the burst such 
as circumburst density and our viewing angle. 

Key words: 


1 INTRODUCTION 


In the las t quarter of 2015 the Advanced LIGO observatory 
fALIGO: [Harrv et al.ll2010ll is expected to begin its first sci¬ 
ence run, spending three months searching for gravitational 
waves (G W). In 2016 it will b e joined by Advanced Virgo 


(AVirgo; lAvernese_^^_^ 2015|), and perform a six -month 


long run I LIGO Scientific Collaboration et al.ll2013l ). Elec¬ 
tromagnetic (EM) follow-up of GW triggers is essential to 
maximise the scientific information gleaned from any such 
detection, as well as providing a useful confirmation of the 
GW signal. To this end, much effort is currently being ex¬ 
pended both in p redicting what EM signal will accompany 
a GW event (e.g. iMetzger Bergerj 12012 : lAndersson et al.l 
12013). and in planning observational foll ow-up strategies 
(e.g. iNissanke. Kasliwal fc Georgievall2013l ). Of course, this 
depends on what gives rise to the GW signal in the 
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ALIGO/Virgo bandpass. While various phenomena (e.g. 
core-collapse supernovae, magnetar flares) are expected to 
give rise to GW signals, the most promising phenomenon is 
believed to be the merger of two stellar mass compact ob¬ 
jects such as neutron stars (NS) or black holes. The final 
stages of such a merger are expected to emit a GW ‘chirp’ 
that could be detected by ALIGO/AVirgo, if the event is 
near enough to Earth. The predicted rates of such mergers 
that will be detectable wi th ALIGO/AVirgo is highly uncer¬ 
tain, with estimates (e.g. lAbadie et al.ll2010l : iGoward et al.l 
l2012t) covering up to 4 orders of magnitude! Electromagnet- 
ically, su ch a merger is accompa nied by a range of emis¬ 
sion (see iMetzger fc Berg^ |2012| ) including that of a short 
gamma-ray burst (sGRB). These are short (^2 s) but ex¬ 
tremely intense (L 10®^ erg s“^) flashes of high en¬ 
ergy radiation, followed by a panchromat ic ‘afterglow’ from 
the shocked circumstellar material; see iBergeil ( 20141 ) for 
a review. The marriage of EM and GW signals from an 
sGRB is extremely valuable as the mass, distance and bi- 
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nary inclination information potentially available from the 
GW signal complements the luminosity, redshift and dura¬ 
tion measurements available from EM observations to con¬ 
strain accurately the energetics and hence physics of the 
merger process. These measurements also have the poten¬ 
tial to constrain cosmology, through combining G W dis¬ 
tance information with sp ectroscopic redshifts (e.g. ISchnt j 
ll986l : lNissanke et al.ll20fg l. Additionally, a precise EM loca¬ 
tion should permit the identification of the host galaxy and 
source environment. 

In this paper we consider ho w best to respond to GW 
triggers with the Swift satellite (iGehrels et al.l l2004 l . We 
are only considering the case that that the GW event is 
accompanied by an sGRB, and we assume that the sGRB 
arose from an NS-NS merger, i.e. we do not consider NS-BH 
mergers. This is primarily because the i nput to our work 
are the GW simulations of ISinger et al.l ll2014l '). who simu¬ 
lated the GW signal from NS-NS systems. The main conse¬ 
quence of this is that we may be underestimating the typical 
distance of the GW triggers, since NS-BH binaries can be 
dete cted out to greater distance than the NS- NS binaries 
(e.g. ILIGO Scientific Collaboration et al.il201^ . Swift was 
designed to study GRBs and, at the time of writing, the on - 
board Burst Alert Telescope IB AT: iBarthelmv et al.lf2005l l 
has detected nearly 1,000, including ~85 sGRB£], and lo¬ 
calised them to within 3' at trigger time. Of those bursts 
which were promptly obser ved by the narrow f ield, 0.3-10 
keV X-ray telescope IXRT: [Burrows et al.|[2005h almost all 
(94% of all GRBs, 80% of short GRBs) have been detected 
and arou nd 40% were detect ed by the UV/optical telescope 
fUVOT: iRoming et al.lBoOSll . Swift-HAT has a peak sensi¬ 
tivity at a significantly lower energy than GGRO-BATSE 
and Eermi-GBM and -LAT, which disfavours the detection 
of sGRBs (the fraction of Swift-triggered GRBs that are 
short is much lower than for these other missions). However, 
very nearby events, such as expected as GW detections, are 
likely to be bright in gamma-rays and largely unaffected by 
this bias. A more significant issue is the BAT’s constrained 
field of view (2 sr, i.e. ~ | of the sky) which means that it is 
rather unlikely that an sGRB which triggers ALIGO/AVirgo 
will also trigger the BAT (this is without considering the 
question of whether typical mergers are likely to produce 
bright sGRBs, which we discuss in Section[7]). Therefore, in 
order to localise the sGRB with Swift we will need to rely on 
follow-up with the narrow-field instruments. In this work we 
focus on the XRT, since it has detected almost every GRB 
it has looked for; and, as we will demonstrate in Section [2. 11 
the rate of contaminating’ X-ray transients unrelated to the 
GW trigger is much lower than at optical wavelengths. 

While the BAT GRB localisation fits easily inside the 
XRT field of view (radius 12.3'); for GW events the error re- 
gion will likely cover many hundreds o f square degrees (e.g . 
Nissanke. Kasliwal fc Georgieval l2013l : ISinger et al.l l201^ : 
Berry et al.l I2OI5I 1. Swift not only has the ability to slew 
rapidly, but it can also auton omously perform multiple 
overlapping ‘tiled’ observations llEvans et al.l l2oi5h , allow¬ 
ing it to cover large error regions, although this has not 
yet been attempted for regions as large as expected from 
ALIGO/AVirgo. Swift has previously observed the fields of 


^ Defined here naively as simply those with Ego < 2 s. 


two GW triggers llEvans et al.ll201^ ’) , but these observations 
were extremely limited, covering less than on e square degree 
between the two triggers. iKanner et al.l ll2012l'l explored ways 
of observing the large GW error regions with XRT, and sug¬ 
gested that, initially, XRT should observe a subset of nearby 
galaxies in the GW error region. When the GW horizon 
distance increases, they suggest observing a series of tiled 
pointings, of 100 s each. Despite this work, it is not obvious 
what the optimal means of observing the ALIGO/Virgo er¬ 
ror region with the Swift-XiKT is, both in terms of where to 
look, and how long to look for. In this paper we use end-to- 
end simulations to explore the optimal approach to following 
up GW triggers with Swift-XiRT, under the assumption that 
the underlying event is a NS-NS merger accompanied by an 
sGRB. 


2 HOW TO IDENTIFY THE X-RAY 
AFTERGLOW 

As noted by authors such as iMetzger fc Berger! 1 I 2 OI 2 I L when 
responding to triggers such as those from ALIGO/AVirgo, 
we do not simply need to detect the EM counterpart, but 
also to id entify it as su c h, fro m among the unrelated sources 
detected. lEvans et al.l (l2015h proposed two ways of distin¬ 
guishing an X-ray GRB afterglow from unrelated, uncata¬ 
logued sources using Swift: brightness or fading behaviour. 
If a source is bright enough that it should have been cata¬ 
logued, but is not, then it is by definition some form of tran¬ 
sient phenomenon, and could be the afterglow. Alternatively, 
GRB afterglows are known to fade on short timescales, 
therefore an uncatalogued and fading object is also a good 
candidate afterglow, even if it is not above the catalogue 
limit. However, GRBs are not the only transient phenom¬ 
ena in the X-ray sky, and when observing the large GW 
error regions we need to consider the possibility of coinci¬ 
dentally detecting a transient, unrelated to the GW signal, 
which we now do. 

2.1 False alarm rates 

To investigate the rate of bright X-ray transients, we se¬ 
lected all GRB datase ts (apart from stac ked images) from 
the ISXPS catalogue ilEvans et al.ll20i3 L which gives po¬ 
sition and flux (among other things) for all point sources 
observed by XRT up to 2012 October. We limited ourselves 
to the GRB fields as, in these cases, which X-ray source was 
the target of the observation is known, so cannot be mistaken 
for a serendipitous or unrelated object. Within each of these 
fields, we searched for X-ray sources which were not the 
GRB, were not catalogued X-ray emitterfl had a detection 
flag of Good, Reasonable or Poor, and were brighter than 
a 3-g upper limit from the Rosat All Sky Survey (RASS; 
IVoees et al.lll999h . calculated at the location of the source 
(see Section [3] for details of the RASS upper limit determi¬ 
nation). That is, we selected all objects which would trigger 
our ‘bright source’ criterion but lie coincidentally in the XRT 

^ Unless the catalogue entry was from one of the earlier XRT 
catalogues, i.e. detected in the same dataset as in the ISXPS 
catalogue. 


© 2015 RAS, MNRAS OOO.ITIITSl 


















































Swift follow-up of ALIGO/AVirgo triggers 3 


field of view. We found 32 such objects, in a total of 52.6 Ms 
of observing time, covering ~ 95 square degrees on the sky. 

To infer the rate of these transients requires knowl¬ 
edge of their duration. If they are all short-lived (com¬ 
pared to the typical duration of an XRT exposure, ~ 500 
s) then the rate at which X-ray transients ‘turn-on’ is 1 
per 1.64 Ms per XRT field of view (0.12 square degrees). 
The observing strategies explored later in this paper cover 
XRT exposure times of 0.006 -1 Ms per GW trigger; with 
those defined as ‘best’ (see Section [5.1115.21) requiring 0.03- 
0.18 Ms per trigger. So, for each GW trigger we follow 
with XRT, this predicts a ~ 2-10% chance of detecting a 
bright X-ray transient unrelated to the GW trigger. This is 
much lower than the expected rates of optical transients, 
which coul d be as high as hundreds of events per GW 
trigger ( e.g. Nissanke, Kasliwal fc Geor gievall201^ . although 
see, e.g.. lGowperthwaite fc Bergerll2015l for wavs to mitigate 
this). It is also an upper limit, since the ISXPS observations 
from which it was derived are of sufficient exposure to detect 
all sources above the RASS limit, whereas, for most strate¬ 
gies explored in Section [3] the exposures are much shorter 
than this, so most of the 32 bright sources in ISXPS would 
not be detected. 

On the other hand, if the transients last long enough to 
all be ‘active’ at the same time G 5-10 years), then we have 
determined that transients exist on the sky with a density 
of 0.3 per square degree. The strategies we will explore in 
this work cover 6-216 square degrees, thus we would predict 
~ 2-65 transients unrelated to the GW event to be detected 
by our observations. The reality is likely somewhere between 
these two values. The Swift data do not allow us to place 
strong constraints on the duration of the 32 transient events 
as, in all but 9 cases, the transient was either already above 
the RASS limit when XRT first observed it, or still above the 
limit when XRT stopped observing it. For the 9 cases where 
both the start and end of the transient phase were seen 
by XRT, the median duration is ~2xl0® s. This is short 
compared to the time interval over which the 32 transients 
were detected (250 Ms), so we suggest that the false positive 
rate is likely to be closer to the first value (one per XRT field 
of view per 1.64 Ms) than the second one (0.3 per square 
deg ree). _ 

iKanner et al.l ll2013l i have also investigated the expected 
rate of X-ray transients, based on XMM-Newton Slew Sur¬ 
vey data. It is difficult to directly compare their results with 
ours, since they used different selection criteria (using only 
the 0.2-2 keV band, selecting sources with a flux 10 x the 
RASS 2-a upper limit, and only those associated with known 
galaxies) which will detect fewer transients than our ap¬ 
proach. They also only gave the spatial density of transients, 
with no temporal information. However, they predict that 
the probability of XRT detecting a transient unrelated to 
a GW trigger in observations covering 35 deg^ and of 100 
s each, is ^1%. If we assume one transient per XRT field 
of view, per 1.64 Ms, we predict 1.8%; these numbers are 
reasonably compatible, given the above differences in mea¬ 
surements, and supports our belief in the lower of the two 
false positive rates we have deduced. 

To determine the rate at which fading X-ray sources 
unrelated to the GW event will be detected is more com¬ 
plicated, since it depends heavily on the observing strategy 
employed. However, we can gain an estimate of the rate from 


I Evans et al.l (l2015l i . In response to 20 neutrino triggers from 
IceCube, they detected a total of 72 sources in their initial 
XRT observations, of which only one initially appeared to be 
fading (further observations showed that is was not related 
to the neutrino trigger). Therefore we can place a crude es¬ 
timate that 1/72=1.4% of objects detected by XRT but 
not related to the GW trigger will show fading behaviour. 

These considerations demonstrate that, if a candidate 
counterpart to the GW event is detected in X-rays, the prob¬ 
ability of it being unrelated to the GW trigger is relatively 
small. Additionally, since Swift is in space with a 96-min or¬ 
bit, and only a 46°-radius exclusion zone around the Sun and 
22'’ zone around the Moon, it can observe a much larger frac¬ 
tion of the sky (and hence GW error region) than a ground- 
based facility. This shows that Swift-XST is a potentially 
valuable resource in searching for the EM counterpart to a 
GW trigger. 


2.2 Feasibility and strategy 


The size of the GW error regions (hundreds of square de¬ 
grees) makes observing them a challenge for the narrow-field 
XRT on Swift', therefore it is necessary to observe a large 
number of ‘fields’: that is, distinct (but possibly overlapping) 
locations on the sky. To determine the optimal strategy for 
observing these fields there are two chief questions which 
need investigating. The first is whether to select XRT fields 
to cover the raw error region produced by ALIGO, or to 
target specific locations within it. With its modest field of 
view, XRT cannot in most cases, expect to cover the entire 
GW error region within a reasonable time after the trigger. 
The simplest approach is, if we can observe N XRT fields, to 
select the N fields which enclose t he most probabilit y (this 
is the concept of ‘searched area’ in ISinger et al.l[2014l i. How- 
ever, sGRBs are expected to exist in or near to ga laxies (e.g. 
iFong. Berger fc Fo:j|20l3 : iTunnicliffe et al.ll2014h . therefore 
it may actually be better to first convolve the ALIGO prob¬ 
ability map with an appropriate galaxy catalogue, and ob¬ 
serve the most probable N fields after this convolution: this 
is the approach which w as followed in earlie r Swift follow-up 
of LIGO -Virgo triggers d Evans et al.l 201^, and was advo - 
cated bv lKanner et al.l (|2012li and iNuttall fc Sutt^ (l2010l i 
for nearby GW events. This approach has the advantage of 
dramatically reducing the area that has to be covered, and 
therefore reducing the average time from the GRB to when 
XRT observes it. On the other hand, galaxy catalogues are 
not complete, and therefore some fraction of the time, this 
approach will mean discarding the XRT field which actually 
contains the GRB. 

The second question is for how long each field is ob¬ 
served. Longer observation^ allow fainter sources to be de¬ 
tected, however they also increase the time it takes to reach 
the field which contains the GRB. Since GRB X-ray af¬ 
terglows fade, this increased delay between the GRB and 
the observation may make the GRB undetectable. Also, the 
longer the exposure, the more unrelated sources we expect 
to detect, making identifying the GRB a greater challenge (if 
we have to re-observe the sources to look for fading, then the 


® Here, and throughout, by ‘observation’ we mean a single expo¬ 
sure of a single field. 
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required observing time scales with the number of unrelated 
uncatalogued sources). 

2.2.1 Galaxy Catalogues 

To investigate the first question - whether or not it is bet¬ 
ter to preferentially observe fields containing known galax¬ 
ies - we require a galaxy catalogue to use in our simula¬ 
tions. In this work w e used the Gravitational Wav e Galaxy 
Catalogue (GWGC; IWhite. Daw fc DhillonI I 2 OIII I : a com¬ 
pilation of various pre-existing catalogues. This gives the 
position, distance, size and orientation of galaxies out to 
100 Mpc. The A LIGO probability maps (Section El are in 
HEALPix format llGorski et al.llioosl l. therefore we created a 
HEALPix-format mask, based on the GWGC. Since NS are 
often given a ‘kick’ at the moment of their birth, and in¬ 
deed sG RBs are often found signif i cantly offset from their 
host fe.g. lFong. Berger &: Foxir2010l : IChurch et aDl201ll ~l. we 
added a 100 kpc halo to each galaxy further than 5 Mpc (for 
galaxies closer than this, the angular projection of a 100 kpc 
halo covers an unreasonably large fraction of the sky, but the 
fraction of binary NS mergers which occur within this dis¬ 
tance is negligible); all pixels in the healpix mask that lay 
within a galaxy or its halo were set to 1, the remaining pixels 
were set to 0. 

The GWGC catalogue is not complete, however an es¬ 
timate of its com pleteness (based on R-b a nd lu minosity) is 
given in fig. 5 of IWhite. Daw fc DhillonI ll201ll '). Based on 
this figure we have tabulated the approximate completeness 
of the catalogue as a function of distance in Table [T] We in¬ 
terpolated within this table to determine the completeness 
of the GWGC at various distances throughout this paper. 
This completeness estimate is based on B magnitude, which 
is likely t o be a reasonable indicator of sGRB rate since, as 
noted bv iFong et al.l (l2013l i. sGRB “demographics are con¬ 
sistent with a short GRB rate driven by both stellar mass 
and star formation.” Therefore, the published GWGC com¬ 
pleteness will give a useful indication of the value or oth¬ 
erwise of using galaxy catalogues in the follow-up of GW 
triggers. 

A drawback to the GWGG is that it only extends to 100 
Mpc, whereas in the 2016 science run, when both ALIGO 
and AVirgo are activ e, we expect to dete ct GW events from 
beyond this horizon llSinger et al.l[20l3l . While alternative 
catalogues exist which do extend beyond 100 Mpc, e.g. the 
2MASS Photometric Redshift Catalog (2MPZ: lBilicki et al.l 
|2014^ - which contains photometric redshifts for galaxies 
detected in the 2MASS infrared survey, with a median of 
2 = 0.08 (~ 300 Mpc) - the completeness of these has not 
been quantified, and the redshifts are photometric, so we 
have not used these in our simulations. 


3 SIMULATION PROCEDURE 

In order to simulate the Swift-yiKT follow-up of sGRBs 
which trigger ALIGO/AVirgo, and explore the impact of dif- 
ferent observing st rategies, we used the GW simulations of 
ISinger et al.l ll2014l 'l. These authors simulated the GW emis¬ 
sion from a large number of NS-NS mergers distributed ho¬ 
mogeneously in space, produced the ALIGO/AVirgo signals 
which would be detected from these and applied the data 


Table 1. Completene ss of the GWGC, based on fig. 5 of 
IWhite. Daw fc DhillonI 1 I 2 OI jl . 


Distance 

(Mpc) 

Completeness 

(%) 

<40 

100 

50 

70 

60 

65 

70 

65 

80 

60 

90 

58 

100 

55 

>100 

Oa 


“ The GWGC only includes galaxies within 100 Mpc, hence the 
sudden cut-off. 

analysis routines to those signals. This was done for the ex¬ 
pected ALIGO configuration in 2015, and for the ALIGO- 
AVirgo configuration (and duty cycles) expected in 2016. 
For the mergers which were detected by this process, the 
probability map and the input binary system parameters, 
are available onlin^. We are only interested in the mergers 
which are accompanied by sGRBs, for which the jet of out¬ 
flowing material is directed towards the Eart h. The opening 
angle of sGRB jets is not well known (e.g. iMargutti et al.l 
I 2 OI 2 I: iFong et aLllioiii . so in order to allow a good num¬ 
ber of simulations, we took a va lue at the wide end of those 
predicted, and selected from the lSinger et al.l ll2014ll simula¬ 
tions all those with a binary inclinatior[3 i < 30°, resulting in 
~ 200 simulated objects each yea r (corresponding to ~ 40% 
of the mergers ‘detected’ in the ISinger et al.l |2014| simula¬ 
tions). [singeT^eU^ (l2014l l provided two probability maps - 
all-sky maps where each pixel’s value is the probability that 
the GW event occurred in that pixel - one map was derived 
rapidly using an algorithm called bayestar, the second cre¬ 
ated using a full Markov Chain Monte Carlo (MCMC) anal¬ 
ysis, whic h takes much longe r to run (typically ~ 109 hr ac¬ 
cording to lBerrv et al.ll2015l L We used the bayestar maps, 
since these are available for all of the triggers on the LIGO 
website (the MCMC map is not), and because, due to their 
rapid creation and their similarity to the full MCMC results 
(for 2015, at least) they are likely to be the maps that will 
be used for real XRT follow-up of a GW trigger. 

As noted above, these simulations were distributed ho¬ 
mogeneously in space, whereas NS-NS mergers are expected 
to occur in or near to galaxies, so we needed to modify the 
spatial distribution of the simulated GW events to reflect 
this expectation. This can be done relatively simply, because 
the shape of the GW error region is determined not by the 
celestial location of the merger, but by the geocentric coor¬ 
dinates. This means that we can shift the simulated mergers 
and probability map in RA without invalidating them, pro¬ 
vided we also also adjust the trigger time such that the GRB 


^ http://www.ligo.org/scientists/first2years/ 

® While binary inclination angle is normally given in the range 
0°—90'’, the GW simulations differentiate between clockwise and 
anticlockwise rotations, allowing i to be in the range 0-180; there¬ 
fore we formally selected objects with i < 30'’ and i > 150'’. 
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is at the same geocentric position as in the original simula¬ 
tions. 

Therefore, for each GW simulation we determined the 
completeness of the GWGC at the distance of the GRB (Ta¬ 
ble [T|). For objects less than 100 Mpc away (i.e. the distance 
range covered by the GWGC), a random number was gener¬ 
ated to determine whether this object should lie in a GWGC 
galaxjQ. If it should, a galaxy was selected at randoirQ from 
those within the GWGC which had a distance consistent 
with that of the simulated merger, and a declination within 
3° of the original LIGO simulatior[3. A position was ran¬ 
domly selected within this galaxy, with uniform probability 
through the galaxy, and then a kick of random magnitude 
(up to 100 kpc) and direction was applied to it. The GW 
probability map was then rotated on the celestial sphere 
such that the GRB lay at this position. For mergers where 
the random number determined that it was not inside a 
GWGC galaxy, or mergers at a distance greater than 100 
Mpc (i.e. not covered by the catalogue) no shift was applied. 

For each merger, regardless of whether we repositioned 
it, the day of the trigger was selected at random. The time of 
day at which the merger would lie at the geocentric location 
used for the GW simulations on this day was determined: 
this was set to be the GRB trigger time. For the 2015 sim¬ 
ulations, the times were all in a 90 day window starting on 
2015 September 1 (i.e. approximately the time expected to 
be covered by the science run); for 2016 we used a 180 day 
window starting in mid June. The only impact this has on 
the simulations is to set the positions of the Sun and Moon, 
hence the regions of sky which Swift cannot observe. 

As well as the GRB location, distance and error re¬ 
gion just described, the following two parameters were also 
needed for each simulation: 

Time delay, Td - the time interval in hours between 
the GW trigger, and the commencement of Swift observa¬ 
tions. Since the initial GW probabili ty map will be pr o- 
duced within minutes of the trigger JSinger et al.l bOlJl 'Fl. 
we drew Td at random from the distribution of delay times 
found in the Swift follow-up of IceCube neutrino triggers 
llEvans et al.ll20i5| j. which is a reasonably analogous ToO 
programme. This is well modelled as a Lorentz distribution: 


P{Td) = 


Ni 


( 

Td-Cl 

V 

Wi 


2 


( 1 ) 


where P{Td) is the probability of the delay being Td, Ni = 
1103, Cl = -1.95 hr and Wi = -8.97 x 10“^ hr. 

X-ray afterglow template. For each simulation 
we selected at random one of the sGRBs listed in 


® i.e. the random number, R, had a value in the range 0 < i? < 1. 
If the GWGC was, e.g. 70% complete at the merger distance, then 
the merger was considered to be in a GWGC galaxy if i? < 0.7. 

^ A better approach would be to weig ht the galaxie s by lu minos- 
ity, however since we cannot move the lSinger et al] ll2014h simu¬ 
lations in RA this is not possible. 

® It is not practical to require the shift to be purely in RA; how¬ 
ever a 3° declination shift is so small that the GW probability 
maps are not invalidated by this shift. 

^ There may be an initial process of human screening before the 
map is distributed, but this is still being debated, so we have 
ignored it. 


iRowlinson et al.l (|2013lf^. for which a valid light curve fit 
exists in the live XRT GRB catalogu^Hl ll Evans et al.ll200^ i. 
We then took the automatic multiply-broken power-law fit 
to XRT light curve from this catalogue, and ‘shifted’ it to 
the redshift of the GW simulation by applying luminosity 
distance and time dilation corrections. For those sGRBs 
where the re dshift was not known, we assumed the mean 
value of 0.72 llRowIinson et allbOHjl^ . We did not apply a 
k-correction, to keep the simulations as simple as possible, 
having first verified that this made little difference: even in 
the most extreme case (taking a short GRB from « = 2.6 to 
simulate a NS-NS merger at 9 Mpc) the k-correction would 
make ^ 10% difference to the light curve normalisation. We 
can only use as templates those sGRBs for which we have a 
valid light curve fit, which introduces a bi as: we have light 
curve fits for 29 of the sGRBs listed in IRowlinson et al.l 
(l2013li (of the remaining 14, 6 were not detected by XRT, 
and for the others we have only a single light curve bin, or 
a bin and an upper limit, so no fit is available). Thus, we 
are really simulating the brightest (in terms of flux, not lu¬ 
minosity) 67% of sGRBs observed by XRT. As will be seen 
in Section [T] this bias is not important when compared to 
other considerations. 

In order to marginalise better over these different pa¬ 
rameters, we performed this process (from deciding whether 
to shift the GW simulation onwards) ten times for each GW 
simulation, giving 2,350 simulations for 2015, and 1,930 for 
2016. 


3.1 Defining the XRT follow-up 

At the present time it is not known how many fields it is 
practical to observe with Swift immediately after a GW 
triggeS nor do we know how many we would choose to 
observe, since a point of diminishing returns is likely to exist 
(i.e. where observing extra fields barely increases the proba¬ 
bility of detecting the GRB). The number of fields observed 
is therefore also a parameter we wish to explore. Initially, 
we created up to 2,400 XRT fields per GW simulation - 
more than we expect to observe in reality - by splitting the 
sky into a pattern of overlapping XRT fields, integrating 
the probability in the GW map within each field, and se¬ 
lected those with the highest probability (up to 2,400 fields, 
or where the probability of the event lying in a field fell 
below 10“®). We did this in two ways: (1) integrating the 
probability from the original GW probability map, and (2) 
by integrating the probability from the GW map convolved 
with the GWGC mask fSection 12. 2. 111 . 

In order to realistically simulate XRT observations, the 
following information was produced for each field, in addi¬ 
tion to the probability of it containing the GRB. 

(i) The mean X-ray background level. This was 


i.e. with a BAT Tqq < 2 s, and to which Swift slewed promptly 
upon detection. 

htt p: / / WWW. swi f t. ac. u k /xrt _Iive_cat 

iD’Avanzo et ah l laoTi) used a slightly different criteria for se¬ 
lecting sGRBs, and found a mean redshift of 0.85. 

Investigative work is under way. 
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Figure 1. The impact of targeting galaxies. When 2,400 XRT 
fields are created to cover the most probable region of the ALIGO 
error region, the solid line shows the cumulative distribution of 
which field (ordered by probability) the GRB actually lies in. The 
dashed line shows the same, if the ALIGO error region is first con¬ 
volved with the GWGG (Section 12. 2. 111 . Top: The 2015 ALIGO 
configuration (based on the BAYESTAR probability map). Using 
the GWGG fewer fields are needed to reach a given probability 
of observing the GRB. Bottom: The 2016 ALIGO-AVirgo con¬ 
figuration: the BAYESTAR-based results are in black, those derived 
from the full MGMG error region are in red. Using the GWGG re¬ 
duces the number of GRBs observed, because most of them occur 
beyond the 100 Mpc limit of GWGG catalogue. 


iKraft. Burrows fc Nous^ lll99ll ). This was converted to an 
XRT count-rate using pimms as above. 

(iv) Details of uncatalogued X-ray sources coinci¬ 
dentally in the field. We derived the se using the log A^-log 
S distribution of iMateos et al.l (l2008l ). which gives the sky 
density of extra-galactic X-ray sources as a function of flux. 
From this we determined the expected number of sources 
unrelated to the GW trigger (Na) in the field of view, with 
a flux of at least 10“^® erg cm“^ s“^ (this limit was chosen 
as it is well below the detection limit of Swift-XKT in the 
exposure durations we are going to explore). Then, drew the 
number of sources (A^s,r) in this simulated field at random 
from a Poisson distribution with a mean of Ns. If this num¬ 
ber was greater than the number of sources already identified 
in step (ii) then we assigned extra sources to the field, until 
it contained (W.r) sources . The fluxes of these sources were 
drawn at random from the lMateos et all ll2008l ) log A(-log S 
distribution, up to a maximum flux equivalent to the RASS 
upper limit for the field. We converted from flux to an XRT 


count-rate using the AGN spectrum defined in step (ii) 


(v) The visibility windows of the field for Swift- 
XRT, for an interval of 15 days starting at the trigger. 


Having thus generated 2,400 fields, the most probable 
N can be selected as inputs to our simulated observations, 
where N is one of the parameters the simulations are de¬ 
signed to explore. 


3.2 Simulating the observations 

Having defined the fields for XRT follow-up, we simulated 
the actual XRT observations for different observing strate¬ 
gies, to evaluate the impact of those strategies on how many 
sGRBs which trigger ALIGO/AVirgo we can identify with 
Swift-NKT. The strategies were based on two parameters: 
the number of XRT fields {N), and the exposure per obser¬ 
vation (e). We explored values of A" = 50,100, 200, 600,1200 
and 1800, and e = 10, 50,100 and 500 s. This was done both 
for the fields based on the raw GW probability map and 
those based on the convolution with the GWGG catalogue. 

We simulated 15 days of observations, recording after 
each field was observed whether the GRB had been identi¬ 
fied, and also how many uncatalogued sources which were 
not the GRB had been identified. Full details of how the 
fields to be observed were selected are given in Appendix 1X1 


drawn at random from th e distribution of s uch levels seen 
in the ISXPS catalogu43 dEvans et al.ll20l3 ). 

(ii) Details of any catalogued X-ray sources in the 
field. These were taken from the RASS. The expected 0.3- 
10 keV XRT count rate was derived from the Rosat count 
rate using pimms and assuming a typical AGN spectrum: an 
absorbed power-law with Ah = 3 x 10^*’ cm“^, F = 1.7. ) 

(iii) The 3-a RASS upper limit in the field. 
This was derived using the RASS images, background 
maps and exposure mapf^. and the Bayesian method of 


1'* We used the 0.3-10 keV background count-rate (‘FIELDBGO’) 
from all ISXPS fields with a ‘FIELDFLAG’ of zero (i.e. they 
contained no artefacts or extended sources). 

Obtained from ftp://ftp.xray.mpe.mpg.de/ftp/rosat/archive 


4 RESULT 1: SHOULD WE TARGET 
GALAXIES? 

4.1 2015 

To investigate the effect of convolving the GW error region 
with the GWGG (Section 12. 2. Ill , we recorded which field 
(ordered by decreasing probability) the GRB actually lay 
in. Not all GRBs lay within the 2,400 most probable fields, 
regardless of whether galaxy targeting was employed. For 
the 2015 ALIGO configuration, the number of bursts which 
did lie in these regions was barely affected by use of the 
GWGG: 1,291 (/2,350) GRBs lay inside the most-probable 
2,400 fields whether or not the GWGG was used; 425 GRB 
only lay in these fields if the GWGG was used, and 394 were 
only present if the GWGG was not used. However, use of 
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Figure 3. The fraction of GRBs identified as a function of how many observations of all fields have been performed. This is for the 2015 
configuration simulations, with galaxy targeting. The 4 panels are for the 4 values of e (exposure per observation) we simulated, and the 
different colours refer to different values of N (the number of fields). 


Solid: 2015. Dashed: 2016 



Figure 2. The distribution of the simulated NS-NS mergers de¬ 
tected by the simulated 2015 ALIGO configuration (solid black) 
and 2016 ALIGO-hAVirgo configuration (dashed red). Both plots 
are normalised to 1. 


the GWGC enables the GRB to be observed quicker. Fig.[T] 
shows the cumulative distribution of which field contained 
the GRB, for both approaches. Gonvolving the ALIGO er¬ 
ror with the GWGG means that we need observe far fewer 
fields to image the GRB. For example, to contain 50% of the 
ALIGO triggers if we use GWGG convolution, we need ob¬ 


serve ‘only’ 500 fields; without the convolution 1,030 fields 
are needed. 

In summary: for 2015 it is best to target galaxies pro¬ 
vided the GW events arise from NS-NS binaries. For NS-BH 
binaries, the horizon distance of aLIGO is greater, and the 
situation is more akin to the 2016 NS-NS case discussed be¬ 
low. 


4.2 2016 

For 2016 the situation is very different. Unlike 2015, most 
of the NS-NS mergers detected were more than 100 Mpc 
away (see Fig. 0 and so do not reside in GWGC galaxies. 
In consequence, targeting those galaxies reduces the num¬ 
ber of NS-NS mergers that lie in XRT fields. Thus for the 
2016 configuration, targeting GWGC galaxies is detrimen¬ 
tal to our chances of detecting the X-ray counterpart to a 
GW trigger. Without galaxy targeting, the typical number 
of fields required to image the merger location in 2016 is 
similar to that in 2015. That is, the 2016 configuration in¬ 
creases the horizon distance without reducing t he size of the 
error region. However, this is a software issue. ISinger et al.l 
1 I 2 OI 4 II show (their fig. 3) that for the 2016 configuration the 
fast BAYESTAR algorithm which is used to construct the GW 
error region does not perform particularly well, and results 
in a larger-than-necessary error region. The full MCMC ap¬ 
proach typically gives much smaller error regions; however 


© 2015 RAS, MNRAS OOO.ITlfTSl 


































































Evans et al. 


Solid: no targeting. Dashed: galaxy targeting 
Black: 50, Red: 100, Blue: 200, Green: 600, Orange: 1200, Grey: 1800 fields 



Figure 4. The fraction of GRBs identified in the 2015 simula¬ 
tions, after a single observation of all N fields, as a function of 
the exposure per observation. The solid lines show the result when 
the N fields are based on the raw ALIGO probability map, and 
the dashed lines for if the map is first convolved with the GWGG. 


it also takes much longe r to run, (‘hours to days’ accord¬ 
ing to ISinger et al■ll2014^■ w hich makes waiting for it im¬ 
practical. Singer et al.l ( 2014h were optimistic that, by 2016, 
the time delay may have been dramatically reduced. Indeed, 
the latest news (Singer, private communication) is that the 
MCMC process has been made faster, and bayestar has 
been improved and now performs as well as the full MCMC 
for the 2016 confi guration. The r efore, the full MCMC prob¬ 
ability maps from ISinger et aP (l2014h are more representa¬ 
tive of what the rapidly available bayestar maps will ac¬ 
tually be like in 2016, so we repeated our simulations for 
2016, this time using the MCMC probability map instead of 
BAYESTAR. Un fortunately, o n ly 100 of the 193 on-axis merg¬ 
ers for 2016 in ISinger et al.l (|2014| ~) had MCMC probability 
maps, resulting in lower statistics. However, Fig. [T] shows 
that, with MCMC, the situation is dramatically improved. 
It is of course still true that use of the GWGC is detrimental, 
since most mergers in 2016 lie beyond that catalogue’s dis¬ 
tance horizon, however even without galaxy targeting, the 
smaller error regions available in 2016 actually mean that 
it takes fewer fields on average to observe the location of 
the GRB than with galaxy targeting in 2015. As noted ear¬ 
lier (Section l2.2.1ll . other catalogues such as 2MPZ may be 
appropriate for use in 2016. 

In summary: for 2016 the GW regions should not be 
convolved with the GWGG. Other galaxy catalogues may 
be appropriate, provided their completeness out to 200-300 
Mpc can be derived and shown to be substantial. 


5 RESULT 2: OPTIMAL OBSERVING 
STRATEGY 

5.1 2015 


lines) and e (different panels) used in our simulations. Recall 
that we only simulated on-axis GRBs, therefore this is the 
fraction of such GRBs, not the fraction of merger events in 
total. For all of the configurations we tested, almost every 
GRB identifiable is identified in the first observation of its 
field: repeat observations have little value, as far as iden¬ 
tifying the bright X-ray object is concerned. The exposure 
dependence is illustrated more clearly in Fig. [d] where we 
show the fraction of GRBs identified after a single observa¬ 
tion of all N fields, as a function of the exposure time per 
observation. Little exposure dependence is seen, except for 
the combination of large values of e and N, where the frac¬ 
tion identified drops slightly. This is because, in such cases 
it takes so long to reach the field containing the GRB that 
fades below the RASS limit before it is observed. Fig. 3] also 
includes the simulations with and without galaxy targeting, 
clearly demonstrating the value of galaxy targeting in 2015. 

In a very small number of the simulations (~ 0.05%) 
the GRB was detected, but was not bright enough to be 
identified. In this case, a second observation could be per¬ 
formed to identify fading behaviour in the source. However, 
this is a very inefficient approach. As Fig. [5] shows, most 
strategies will result in the detection of new, uncatalogued 
sources which are not the GRB, all of which would have to 
be followed up; and as already noted, the number of cases 
where the GRB is among these sources is extremely small. 

5.1.1 False positives 

We noted earlier (Section 12.111 that the expected rate of X- 
ray transient sources is 1 per 1.64 Ms of exposure time. So 
for the case of 600 fields, observed for 50 s each, there is 
only a 2% chance of serendipitously observing a X-ray tran¬ 
sient unrelated to the GW trigger. There is, however, an ad¬ 
ditional source of potential contamination: due to Poisson 
noise, an uncatalogued source with a flux below the RASS 
limit will sometimes be measured to have a flux more than 
1 -cr above this limit, thus meeting one of our criteria to 
be identified as the likely GRB afterglow. Our simulations 
included Poisson noise, and we recorded how many of the 
unrelated sources detected would be thus classified as a pos¬ 
sible afterglow. Fig. |6] shows the fraction of simulations in 
which this happened for at least one unrelated source. This 
increases sharply for exposures longer than 50 s. This means 
that we cannot interpret the discovery of a bright (>RASS 
limit), uncatalogued source to be an unambiguous detection 
of the sGRB: it could simply be a fainter source subject to 
Poisson errors. Also, there could be multiple ‘bright’ sources 
detected in the follow-up of a single GW trigger. Therefore, 
it will be necessary to re-observe any such sources for ~ 1 
ks each to confirm which (if any) of them was a transient 
event. The maximum number of such false positives detected 
in any of our simulations is 19, so this confirmation obser¬ 
vation requires only a small investment of exposure time, 
compared to that used to detect the sources. 


The previous section demonstrated that, for the 2015 con¬ 
figuration, it is preferable to target galaxies in the ALIGO 
error regions. In Fig. |3]we show the fraction of the galaxy¬ 
targeting simulations for which the GRB was identified, as a 
function of how many observations of the N fields were con¬ 
ducted. This is shown for all values of N (different coloured 


5.2 2016 

We performed two sets of simulations for 2016, one based on 
the BAYESTAR error maps and one based on the full MCMG 
maps. Since Fig. [T] showed that targeting GWGC galax¬ 
ies was detrimental in 2016, we only simulated observations 
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Figure 5. The cumulative distribution of the number of uncata¬ 
logued sources unrelated to the GW trigger detected in the 2015 
simulations (with galaxy targeting), assuming observations of 50 
s (top panel) or 500 s (bottom). If no bright source is detected in 
the observations, all of these uncatalogued sources would have to 
be reobserved to search for signs of fading. 



Figure 6. The fraction of simulations in which a source unrelated 
to the GW trigger was classified as a potential afterglow, because 
its measured flux was above the RASS limit due to Poisson noise. 
The data shown here are for the 2015 simulations, with galaxy 
targeting. 

based on the raw GW error maps. As with the 2015 simula¬ 
tions almost all identifiable GRBs were identified in a single 
observation. Fig. Q shows the fraction of GRBs detected as 
a function of exposure time, for the different simulations, an 
analogue of Fig. |4] Using the bayestar maps we identified 
a lower fraction of GRBs than in the 2015 simulations, al- 


Solid; bayestar. Dashed: MCMC 

Black: 50, Red: 100, Blue: 200, Green: 600, Orange: 1200, Grey: 1800 lieids 



Figure 7. The fraction of GRBs identified in the 2016 simula¬ 
tions, after a single observation of all N fields, as a function of 
the exposure per observation. The solid lines show the result if 
the BAYESTAR error map is used, and the dashed lines are for the 
full MGMG error map. 

though for the N < 200 simulations, the difference is less 
obvious: the fraction identified in 2016 is comparable to the 
2015 follow-up with no galaxy targeting. In contrast, using 
the full MCMG-based error map in 2016 we identified the 
GRB in a larger fraction of simulations than in 2015. 

One curious feature of Fig.[7]is that, for the N — 1,800 
simulations based on the full exposure map, fewer GRBs 
were identified than the N = 1, 200 simulations at e = 500 s. 
This is an artefact of the way the simulations decide in which 
order to observe the XRT fields (Appendixl^, which means 
that while the fields in the N = 1,200 simulation are the 
most probable 1,200 fields in the N = 1, 800 simulation, they 
are not necessarily the first 1,200 to be observed in the latter 
simulations. This, combined with the spatial distribution of 
probability in the 2016 MGMC error maps, which tends to 
contain distinct ‘clumps’, means that on average the field 
which actually contained the GRB was observed later in the 
N = 1,800 simulations than in the N = 1,200 case. As a 
result, the GRB was more likely to have faded to below the 
detection threshold or RASS limit. 


6 IF FERMI TRIGGERS AS WELL 

If Swift-HAT and ALIGO/AVirgo trigger on the same event, 
the discussion in this paper is moot: the BAT error region 
fits comfortably in the XRT (and UVOT) field of view, and 
the satellite will perform the standard automatic follow-up 
process which has resulted in the rapid detection of the 
X-ray afterglow of almost all sGRBs detected by Swift- 
BAT to date. However, as noted earlier, this is unlikely. 
A more probable scenario for the simultaneous GW and 
EM detection of the same event is for Fermi to trigger on 
the GRB. Fermi conta ins two instruments : the Gamma-ray 
Burst Monitor ( GBM: iMeeean et al. l2009l) and Large Area 
Telescope (L AT; I Atwood et al.l |200£ ). GBM is sensitive to 
triggers over the entire sky which is unocculted by the Earth, 
although has poor positional accuracy, with typical uncer¬ 
tainties of many degrees. LAT has a smaller field of view 
(2.4 sr), and if GBM detects a bright GRB Fermi with 
a sufficiently high peak flux, Fermi automatically slews to 
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Black: 2015 (solid=aLIGO error. Dashed=galaxy targeting) 
Red: 2016 {solid=bayestar. Dashed=full MCMC) 



Field containing the GRB (ranked by probability) 

Figure 8. As Fig. ^ but showing the effect of using informa¬ 
tion from a Fermz-GBM trigger simultaneous with the GW trig¬ 
ger. The black lines are for the 2015 ALIGO-detected mergers, 
with the dashed and solid lines showing respectively whether 
the GWGC is or is not used. The red lines show the 2016 
ALIGO/AVirgo objects, with the solid line showing the bayestar 
probability maps and the dashed line the full MCMC map. 


move the GRB into the LAT f ield of view (|Ackermann et ahl 
I2013I: I von Kietilin et al. llHi). LAT positions can be much 
better than those from GBM, down to While it is 

beyond the scope of this work to simulate the probabilities 
of these instruments triggering on GW-detected sGRBs, we 
consider how such triggers will influence Swift observations. 
Note that, due to their proximity, we would expect such 
sGRBs to be bright enough to trigger Fermi (or indeed BAT 
or the Interplanetary Network [IPN]) unless they are ~ 100 
times fainter than the typical Swift-detected sGRBs (see 
Section III). 


6.1 If LAT triggers 

While the LAT error region is highly variable, at the lower 
end they are small enough to fit inside a single XRT field 
of view. For less-well constrained positions, Swift has often 
observed the LAT error regio n in a small number of tiles 
(typically 4, sometimes 7; e.g. lEvan 3 l2013h . This allows us 
to rapidly cover the entire error region, and build up long 
exposures of multiple kiloseconds per tile in much less time 
that it takes to observe each field once if we have only a GW 
trigger. These longer exposures allow easy identification not 
just of a bright source, but also of a faint, but clearly fading 
object. Due to the small area covered, the possible presence 
of unrelated sources is not a signihcant issue, since the bur¬ 
den of following multiple sources to confirm their behaviour 
is much reduced when those sources are close together. A 
drawback to LAT triggers is that nearly all of them are 
found in ground processing, which means that are not de¬ 
tected until many hours after the event. However, a nearby 
burst will likely also be bright, increasing the likelihood of 
an onboard LAT trigger. Further, even if Swift has already 
started tiling the GW error region, should a LAT error re¬ 


gion become available, we could immediately terminate the 
ongoing observations and point Swift at the LAT region. 


6.2 If GBM triggers 

Even if LAT does not detect the GRB, a GBM detec¬ 
tion simultaneous with an ALIGO/AVirgo trigger has the 
potential to significantly reduce the sky area we need to 
search. The median 1-a statistical error radius for a GBM- 
detected sGR.iIR is 7°. A systematic error is also present 
(jConnaughton et al.ll201,^ L which we simplistically treat as 
a Gaussian with cr=3.7°, giving an overall error median er¬ 
ror of 8°. For each simulated GW event in our samples, we 
selected at random the location of the centre of the GBM 
error region, assuming a circular region, with a Gaussian ra¬ 
dius, with (Ttot = 8°, and then created the list of XRT fields 
to observe as in Section^ but also only accepting fields that 
lay within S-ctot of the simulated GBM position. We then 
determined how many fields we would have to observe (to a 
maximum of 2,400) in order to observe the GRB location. 
The result is shown in Fig. [8] which can be compared with 
Fig. [T] This shows, not surprisingly, that if a GBM locali¬ 
sation is available, it can reduce the number of XRT fields 
we need to observe significantly, e.g. in 2015, with galaxy 
targeting, the number of XRT fields needed to observe 50% 
of the GRBs is just 190 (compared to 500 without GBM: 
see Section . 


7 ARE DETECTABLE, NEARBY SGRBS 
LIKELY? 

So far we have focussed on how best to search for an XRT 
counterpart to a GW trigger, but another important con¬ 
sideration is whether a GW trigger from a nearby, on-axis 
sGRB is likely. Various works have attempted to calculate 
the rate of nearby sGRB s, either from the t heoretically ex¬ 
pected merger rate (e.gjAbadle_^£_aD|20]^ or the rate of 
observed sGRBs (e.g. ICoward et al. 20121 ). As noted ear- 
lier, these est i mates have large uncertainties, for example, 
ICoward et al.l (l2012l) give an sGRB r ate d ensity of 8 - 
1100 Gpc"® yr-^ 


Metzger fc Berg^ 1I2OI2I ) fitted the ob¬ 


served redshift distribution of S'uij/t-detected sGRBs with 
known redshift, divided this by the fraction of sGRBs with 
known redshifts, and from that inferred that “a few Fermi 
bursts over the past few years probably occurred within 
the ALIGO/Virgo volume.” However these approaches do 
not consider one important observational constraint: unless 
the nearby sGRBs are systematically less luminous than the 
ones for which we have redshifts, they should appear much 
brighter than those - orders of magnitude brighter. Is there 
evidence for any such bright objects among the sGRBs de¬ 
tected to date? 

To investigate this we want to compare the expected 
brightness of an sGRB inside the ALIGO/Virgo horizon 
with the distribution of observed sGRB brightnesses. We 


Taken from _ the Fermi _ GRB cat alogue , 

iGruber et alJ (120141 ): Ivon Kienlin et ahl |2 m3), 

http://heasarc.gsfc.nasa.gov/W3Browse/fermi/fermigbrst.html, 
selecting objects with Tgo < 2 s and flux_batse_64> 2.37 s~^, 
following I Wanderman Sz PiranI |201^. 
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Figure 9. A comparison of the observed brightness of the Swift-sGKBs without redshift, to what would be expected from an sGRB 
within the GW horizon distance. Pq 4 is the peak 64-ms photon flux. The blue histogram shows the distribution of Pg 4 for the Swift 
short GRBs without redshift. The red histogram shows the distribution of short GRBs with redshift measurements, with the peak flux 
corrected to a distance of 100 (left panel) or 200 (right panel) Mpc. The resultant peak fluxes are ord ers of magnitude larg e r tha n 
anything Swift has observed to date. The solid black line shows the theoretical sGRB luminosity function of lWanderman fc Piran ll201fil . 
converted to P64 a^t the two distances. 


define the brightness as the 64-ms p eak photon finx (P 64 ) 
from the prompt emission, following IWanderman fc PirarJ 
(l2015l) . For the Swift-detected sGRBs with known redshift 
listed in IWanderman &: PirarJ ll2015l ). we derived the 15- 
150 keV P 64 values that would have been measured had 
the GRB been at 2 = 0.025 (~ 100 Mpc) and 2 = 0.05 
(-^200 Mpc); these are shown as red histograms in Fig. [O] 
This sample is likely biased towards intrinsically brighter 
objects (i.e. objects which are easier to detect), so we have 
also convert ed the theoretical Inminosit y function of short 
GRBs from IWanderman fc PirarJ (l2015l ) into the distribn- 
tion of 15-150 keV P 64 at these two redshifts; this is shown 
as the black line in Fig. [O] For comparison, we also show 
the distribution of the observed 15-150 keV P 64 valnes for 
the Swift-BAT detected sGRBs without a known redshift 
(blue histograms). The P 64 distribution of the Swift sGRBs 
with no redshift is clearly not consistent with them being 
identical to the known-redshift sample bnt lying at low red¬ 
shift: in this case the blue histograms would appear sim¬ 
ilar to (or at least share significant overlap with) the red 
ones. It is possible that a small number of the no-redshift 
sample were at ~ 200 Mpc, however they must have been 
of lower intrinsic Inminosity than the with-redshift sample, 
since they just barely overlap with what the latter sample 
would have looked like at 200 Mpc. We certainly cannot rule 
out that some of the no-redshift samp le lay at 100-200 Mpc 
and ca me from the faintest end of the IWanderman fc Pirar] 
(I2OI5I) luminosity function, although this faint end is not 
well constrained. There are, however, two objections to this 
premise. First, the observed P %4 distributions of the with- 
and without-redshift Swift sGRB samples are very similar 
(Fig. 1101) : a Kolmogorov-Smirnov test gives the probability 
that the two samples are drawn from the same distribution 
as 0.43; so if the no-redshift sample does indeed contain 
bursts that are of lower luminosity than the with-redshift 
sample, the redshift distributions of the two samples have 
fortuitously been tuned correctly to give apparently identi¬ 
cal P 64 distributions, which seems a little far-fetched. Sec- 


Black: No redshift. Red: With Redshift 



(ph cm-2 s-') 

Figure 10. The distribution of observed peak 64-ms photon flux 
from the Swift GRBs with redshift (red) and without (black). The 
K-S test gives a 0.43 probability that the two samples come from 
the same parent distribution. 


ondly, if any of the no-redshift sample really were within 
200 Mpc, we would expect the host galaxy to have been 
found in optical follow-up. and hence the GRB to have a 
redshift. iTunnicliffe et all(120141) explored this question and 
found tentative evidence that one or two ‘hostless’ sGRBs 
could have arisen from stars kicked out of nearby (< 200 
Mpc) galaxies but that these would have surprisingly low 
jet energies. For example, they found a galaxy at a distance 
of 80 Mpc, which could be the host of GRB 111020A, but 
this would require the isotropic equivalent for that GRB to 
be just 10"*® erg, several orders of magnitude lower than typi¬ 
cal sGRBs. These considerations suggest that no compelling 
evidence exists for a BAT-detected sGRB within 200 Mpc. 

Swift-BAT detects fewer sGRBs than GBM does or 
BATSE did, so we have repeated the above test, for the 
BATSE and GBM bursts. The results are shown in Fig. [11] 
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Figure 11. As Fig. [9] but for the 50-300 keV band. The blue and grey histograms show the distribution of observed peak 64-ms photon 
flux for the BATSE and GBM sGRBs respectively — these have no measured redshift. 


For these, we derived P&a in the 50-350 keV band used 
by those instruments. We again plotted the expected P 64 
values for the with-redshift Swift sGRBs (corrected to the 
50-350 keV band), had they been at 100 or 200 Mpc (red 
histograms). The BATSE and GBM observed P 64 distribu¬ 
tions are shown in the blue and grey histograms respectively. 
As with the Swift data, we cannot rule out the possibil¬ 
ity that a few of the sGRBs detected by BATSE or GBM 
could have been within 100-200 Mpc from Earth, however, 
again, these must have been significantly fainter (typically 
at least an order of magnitude) than the 5'uii/t-detected 
GRBs with redshift. A substantial fraction of the observed 
GBM and BATSE bursts are consistent with being at 200 
Mp c, provided they com e from the low-luminosity end of 
the IWanderman fc PiranI ll2015ll luminosity function. How¬ 
ever, this cannot be systematically the case, since the GBM 
and BATSE samples must contain the more luminous GRBs 
that Swift has detected as well. 

Two historic events warrant explicit mention: GRBs 
051103 and 070201 were both very bright events, localised 
by the IPN, and with large error regions that partially inter¬ 
sected M81 and M31 respectively. However, LIGO observa¬ 
tions showed that GRB 070201 could not have been a NS-NS 
or NS-BH merger in M31 ll Abbott et al.l 1200^ 1. suggesting 
instead that it was a giant flare from a soft gamma repeater 
(SGR), or a more distant sGRB coincidentally close in pro¬ 
jection to M31. The case for GRB 051103 is less clear, and 
an SGR origin of this event cannot be ruled out. If it were a 
sGRB in M81 then it had a very low energy release (isotropic 
equivalent of ~ 10^ ^ erg) and an extremely faint afterglow 
dHurlev et aLllioi^ l. 

These considerations mean that we cannot conclusively 
rule out that some of the observed GRBs without redshifts 
could have been nearby - particularly those detected by 
BATSE, GBM or IPN, where the error regions are too large 
for host galaxy searches. However, we can state that, were 
this so, these nearby bursts were also underluminous com¬ 
pared to the with-redshift sample of Swift sGRBs. In more 
than 20 years of observations, no sGRB has been detected 
consistent with being less than 200 Mpc away and having a 
prompt luminosity typical of the Swift sample. This means 
that events such as those we have simulated in this work are 


rare. So if ALIGO/Virgo does detect a binary NS merger, 
it is likely to be significantly fainter (1-2 orders of magni¬ 
tude) than in our simulations, and therefore much harder 
to detect. More likely is that a NS-NS merger detected by 
ALIGO/AVirgo will be off-axis, that is, the jet will not be 
pointed towards Earth. For such objects no prompt gamma- 
ray emission is seen from Earth, so such objects are not in 
our observed GRB sample. Afterglow emission is still ex¬ 
pected from some such objects however, and we now briefly 
consider how detectable such emission may be. 

7.1 Off-axis GRBs 

Although there are no X-ray observational constraints on 
the properties of GRBs viewed off-axis (by which we mean 
that the binary inclination i is larger than the jet open- 
ing angle, Oj), there have been t heoretical predictions (e.g. 
Ivan Eerten fc MacFadv^ l201ll : iGhirlanda et alll2015l j. In 
order to give some indication of the potential observabil¬ 
ity of off-axis sGRB afterglows with Swift we have taken a 
highly simplified approach. First, we assumed that the time 
at which the afterglow becomes visible off-axis (toa), for a 
binary at inclination i radians, is the same as the time at 
which a jet-break would occur if the jet opening angle {6j) 
was i'. 



where E 52 is the energy in the jet (Ejet) in units of 10®^ erg, 
no is the circumburst density in units of cm“®, and toa is in 
the rest-frame of the GRB. 

Rather than the generous value of 30° used earlier, we 
bega n by assuming a jet opening angle of 1 0° for all bursts 
(e.g. IPuffell. Quataert fc MacFadven l2015h . To simulate a 
GRB, we assigned it a distance D up to 200 Mpc (with 
probability P{D) oc D^) and a binary inclination i in the 
range 0-90° (with probability P{i) oc sini). We selected 
an observed short GRB X-ray afterglow to use as a tem¬ 
plate, as in Section |3] but modified the light curve slightly 
to include a jet break (the jet break time calculated as in 
equation (2] but replacing i with 6j = 10°), after which the 
light curve decayed with a power-law index of —2.5. For 
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Figure 12. The number of bursts (out of 100,000 simulated) detectable in a 50 s (black) or 500 s (red) observation, as a function of 
how long after the merger the observation begins. The solid lines are bursts viewed ‘on-axis’, that is the 10° jet is oriented towards the 
observer. The dashed lines are the off-axis events. The different panels correspond to different values of jet energy (Fjet) and circumburst 
density (no). At early times the on-axis bursts dominate the visible population, despite comprising only 1.5% of all mergers. At late 
times, as the afterglow becomes visible off-axis, the off-axis population starts to dominate. 


on-axis binaries {i < 10°) the brightness of the GRB af¬ 
terglow at a given time was taken from this light curve. 
For off-axis binaries, we assumed there was no X-ray emis¬ 
sion before toa, after which the brightness followed the on- 
axis (post jet-break) light curve just defined. Comparison 
with the simulated light curves of ivan Eerten &: MacFadvenI 


ll201llf1 shows that this simplification is not unreasonable. 
We then calculated whether the GRB would be detected in 
observations of 50 s or 500 s, as a function of the time since 
the merger event (for these simple simulations we assumed 


http://cosmo.nyu.edu/afterglowlibrary/sgrb2011.html 
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Figure 13. How many on/ofT-axis afterglows are detected as a function of time for jet angles of 20° (top panels) and 30° (bottom panels). 
Rather than reproducing all the panels of Fig. [TJ] we show only the densities which give the most off-axis bursts for each of the energies 
investigated. The solid lines are bursts viewed ‘on-axis’, the dashed lines are the off-axis events. 


the XRT background level was always 10 ^ ct s ^ pixel 
a typical value from lSXPS: lEvans et al.ll2ni3) . 

We simulated 100,000 GRB s for each combination of 
i?52=0.01 and 0.0001 (following Ivan Eerten fc MacFadvenl 
[2^3), and no=of 0.01, 0.1 and 1 cm“® (e.g. iBerged 2014h . 
Note that, since Ga and 6j scale with E^ 2 lno, some of these 
combinations are degenerate. 


In Fig. [12] we show the number of simulations in which 
the GRB was detectable in observations of 50 s or 500 s, as a 
fnnction of the time at which the observation began. The on- 
axis and off-axis GRBs are differentiated. Although there are 
many more off-axis bursts (98.5% of those simulated), they 
are less detectable, since, by the time the emission becomes 
visible to the observer, it has faded significantly - Fig. 1141 
shows the cumulative distribution of binary inclinations for 
the GRBs detected and shows that this still corresponds to a 
fairly restrictive range of viewing angles. Fig. [T^ shows that 
the detectability of the off-axis objects depends strongly on 
the jet energy and the circumburst density; in reality some 
(not currently constrained) distribution of these parameters 
is expected, which means that the ‘true’ shape of Fig. ll2l for 
the short GRB sample is some form of weighted average of 
the panels shown. 

This variation in detectability makes it difficult to de¬ 
fine the optimal time to observe in X-rays in order to de¬ 
tect off-axis emission. For on-axis GRBs, it is clear that the 



Binary inclination (“) 


Figure 14. The cumulative distribution of binary inclinations 
for the mergers which were detected by our simulations, when 
the jet opening angle was 10°. The range of angles which result 
in detectable GRBs depends strongly on the jet energy and cir¬ 
cumburst density. 


earlier the observations, the more GRBs are detected. For 
off-axis objects, observing too early prevents detection (the 
afterglow emission is still beamed away from the observer). 
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as does observing too late (the emission is too faint), but the 
parameters of the window in which the afterglow can be de¬ 
tected vary according to the burst properties. Fig. ll2l shows 
that by ~ 36 hours after the trigger, the number of off-axis 
bursts detectable exceeds the number of on-axis bursts for 
all Ez 2 = 0.0001 simulations. For if E ^2 = 0.01 then on- 
axis bursts dominate until at least 72 hours post merger. 
A sensible Swift response to a GW trigger may therefore 
be to begin observations as soon as possible (in case the 
merger was on-axis), but after ~36 hours, if the afterglow 
has not been detected, the strategy should instead turn to 
looking for off-axis afterglows which are starting to emerge. 
In this case, rather than continuing the XRT observing pro¬ 
gram (which will be moving to lower-probability portions of 
the GW error region) it would be preferable to re-observe 
the high-probability fields already studied. This has an ad¬ 
ditional science bonus, in that if the GRB afterglow is de¬ 
tected in these second observations, but was not present in 
the first, it confirms that the merger was off-axis, and allows 
us to place constraints on i — Oj. One could envisage a third 
set of observations beginning at 72 hours after the trigger 
(i.e. where the off-axis start to dominate over the on-axis 
GRBs if i ?52 = 0.01). However at this time, the number of 
off-axis afterglows detectable if E 52 = 0.0001 is still much 
higher than if E52 = 0.01, therefore unless there are com¬ 
pelling reasons to believe that the true sGRB population is 
dominated by the higher-energy objects, it is better to con¬ 
tinue the observations than the return the the previously 
observed-fields. 

We also see from Fig.[T2]that the choice of exposure time 
is much more important for off-axis GRBs than for on-axis 
GRBs, with the number of off-axis GRBs detected a factor 
of 3-5 higher for 500 s exposures than for the 50 s exposures: 
for on-axis GRBs the difference is at most a factor of 1.5. 
Therefore, if this ‘double-headed’ approach were employed 
it would make sense to use shorter exposures in the first 36 
hours, while searching for on-axis objects, and then switch 
to longer exposures when the search focus shifts to off-axis 
GRBs. 

While the jet opening angle of 10° used above 
may be reasonably representative of typical jets (e.g. 
iDuffelL Quataert fc MacFadvenll2015^ . l ower limits as large 
as 25° have been deduced for sGRBs (iGruoe et al.l 120061 : 
IZhang et al.l l2015l l. We repeated our simulations for jet 
opening angles of 20° and 30°, extending the runs to 28 days 
after the merger, to account for the later rise time of the off- 
axis jets. Fig. 1 131 shows the results for the values of i? 52 /no 
which maximise the number of off-axis bursts detectable. In¬ 
creasing the jet angle had two effects: some off-axis GRBs 
became on-axis GRBs, and the rise time of the off-axis GRBs 
moved to later times. The latter also means that the off-axis 
afterglows are fainter. These factors conspire to make the off- 
axis afterglow less significant as a proportion of afterglows 
detectable with XRT. 

All of the considerations discussed in this section show 
that chasing off-axis afterglows with a narrow-field instru¬ 
ment such as Swift-XKT is extremely challenging. The num¬ 
ber of such afterglows that can be detected, and the time 
at which they are detectable, are both highly variable, de¬ 
pending on the details of the jet and the medium in which 
burst occurs. 


8 DISCUSSION AND CONCLUSIONS 

We have performed simulations of Swift-XRT follow-up of 
GW triggers from the ALIGO/AVirgo observatories, accord¬ 
ing to their expected performance in the 2015 and 2016 sci¬ 
ence runs. In each case, we assumed that ALIGO/AVirgo 
detected a NS-NS merger accompanied by an sGRB with 
an X-ray afterglow similar to those observed by Swift-yiRT 
to date. We have found that XRT can identify the GRB in 
up to ~ 30% of cases, depending on the observational strat¬ 
egy employed. We considered this strategy in terms of two 
questions: which portions of the GW error region to observe, 
and how best to observe them. 

For the former question, it is clear that in 2015 preferen¬ 
tially observing galaxies in the GW error region significantly 
increases the fraction of GRB afterglows (arising from NS- 
NS mergers) that we can identify with XRT. Beyond 2015 
(or when considering NS-BH mergers), as the sensitivity of 
the GW network increases, galaxy catalogues that extend 
to greater distances (200-300 Mpc) while remaining highly 
comple te will be necessary . One possibility is the 2MPZ cat¬ 
alogue (iBilicki et al.ll20i4 1: its advantages are large sky cov¬ 
erage and uniformity, its disadvantages are use of photomet¬ 
ric redshifts and the current lack of detailed completeness 
analysis. 

If Fermi-GBM also triggers on the event detected by 
the GW facilities, we can convolve the GBM and GW error 
regions, which reduces the amount of sky we need to cover, 
typically by a factor of ~ 2. 

Goncerning the optimal strategy, our simulations show 
that the fraction of GW-detected sGRBs we can detect with 
Swift-ySKT does not strongly depend on the amount of time 
for which we observe each field in our follow-up. Not sur¬ 
prisingly, the more fields are observed, the more GRBs are 
identified, but this is a sub-linear relationship: doubling the 
number of fields less than doubles the number of identified 
counterparts. 

As we pointed out in Section [3] because we used the 
Swift X-ray afterglows as templates in our simulations, our 
results are biased towards the brighter sGRBs. Further, 
we have argued (Section [Tj that sGRBs like those seen by 
Swift occur only extremely rarely within the distance that 
ALIGO/AVirgo will probe in the next few years, suggesting 
that if a NS-NS merger is detected by these facilities, the 
accompanying sGRB will either be (much) less luminous 
than those we have simulated, or will have its jet oriented 
away from us. I ndeed, if sGRBs follow a lum inosity function 
such as that of lWanderman fc PiranI ll2015ll . a nearby GRB 
will likely be low luminosity compared to the ones Swift has 
detected to date. These considerations suggest that longer 
XRT exposures are preferable, to increase our sensitivity 
to these fainter signals: in Section 17.11 we showed that the 
number of off-axis objects detectable increases by a factor 
of 3-5 when 500 s exposures are used, compared to 50 s 
exposures. However, fainter afterglows will not be as easy 
to identify from among the uncatalogued unrelated sources 
also detected in longer exposures. In this case, to identify the 
GRB requires a second observation, to search for evidence 
of fading. This observation will have to be performed for 
every uncatalogued source detected, since we cannot know 
a priori which (if any) is the GRB. As the lower panel of 
Fig. [ 5 ] shows, for the longer initial observations (500 s per 
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field), this may be many hundreds of sources. lEvans et al.l 
ll2015l) showed that typically an exposure of 1 ks would be 
enough to identify whether a source detected in these initial 
observations had faded, with S-cr confidence, so a simplistic 
calculatiorF^ based on Fig. [5] suggests that the time needed 
to perform these follow-up observations is similar to the time 
taken to perform the initial observation of the GW error re¬ 
gion. 

Even if we only perform short exposures with Swift, i.e. 
focus purely on the bright but rare events which are de¬ 
tected above the RASS limit, some false positives will still 
be detected simply because of Poisson effects. Therefore, it 
will be necessary to re-observe even the bright, uncatalogued 
sources; however there are many fewer of these than of the 
faint uncatalogued object just discussed so the required ob¬ 
serving time is correspondingly lower. 

If we do perform follow-up observations of any uncat¬ 
alogued sources, and identify fading behaviour, we cannot 
rule out that the object is an unrelated variable object. As 
we noted in Section [2. II roughly ~ 1/72 sources unrelated to 
the GW event are expected to show strong signs of fading, 
therefore if a fading source is found, continued observations 
will be necessary to definitively identify it as the GRB af¬ 
terglow. 

The discussion so far has focussed on the attempt to 
identify the counterpart to the GW trigger based on Swift 
data alone, whereas a large number of facilities spanning 
the EM spectrum have committed to the follow up of GW 
triggers. While X-rays may represent the best wavelength at 
which to rapidly detect an on-axis sGRB afterglow, sGRB 
emission is panchromatic, and for off-axis events the optical 
and infrared emission from a kilo nova represents arguabl y 
the best hopes for a detection fsee lMetzger fc Bergerll2012l) . 
Kilonova emission is powered by the radioac tive decay of r- 
process elements created during the merger (R osswog et al.l 


I 1999 I : iFreiburghaus. Rosswog fc ThielemannI 19991 ). and de¬ 
tection of such emission has been claimed for GRB 130603B 
llTanvir et al.ll2013l : lBerger. Fong fc Chornockll2013l), as well 
as possibly GRB 0 80503 JP erlev et al.ll2009l : Gao et aD2015l ) 
and GRB 060614 (lYang et al.ll2015l) . The knowledge of what 
other observatories have detected - or have not detected - 
will therefore be of significant value in trying to establish 
which, if any, of the X-ray sources detected is the GW coun¬ 
terpart. However, signals such as the kilonova may not rise 
until some time after the trigger in which case these data 
would not be available to us at the time of the Swift obser¬ 
vations, so the ability to identify the counterpart from Swift 
data alone is desirable. 

Due simply to geometric effects, off-axis mergers are 
much more likely than on-axis ones: for a 10° jet, only 1.5% 
of mergers are on-axis. However, a much lower fraction of 
off-axis objects are detectable by XRT, and only within a 
specific window of time after the merger event. The tim¬ 
ing of this window depends strongly the jet opening angle, 
the energy in the jet, and the density of the circumburst 
medium. The distributions of these parameters for sGRBs 
are not well known, therefore it is difficult to identify the 
‘best’ time to observe when searching for an off-axis after- 


Neglecting, for example, the possibility of observing multiple 
sources in a single field of view. 


glow. This in particular highlights the fact that the ideal 
facility would be a wide-field (prefera bly all-sky) X-ray im - 
ager with good position accuracy (e.g. lOsborne et ^l2013l ). 

As Swift-yiRT has a narrow field of view, we require 
many pointings to observe the GW error region, and new 
flight software to support the upload of hundreds of ‘Auto¬ 
matic Targets’ is being developed to support this. We have 
simulated some fairly extreme observing programmes: for 
example, having Swift observe 1,800 fields for just 10 s each; 
whether this is practical is currently unclear and, at the time 
of writing, the maximum density of observations with Swift 
is still under study. However, we do know that the most 
slews Swift has performed in a single day is 167. This cor¬ 
responds to just over 500 s per observation, if we ignore the 
time spent slewing or passing through the South Atlantic 
Anomaly (when XRT cannot collect data). Therefore, our 
500-s per observation simulations give us a reasonable idea 
of what we can expect in reality. In this regime, if we ob¬ 
served for 1 day (i.e. ~ 150 fields) Fig.|3]shows us that, if the 
NS-NS merger were on-axis, we would identify the sGRB in 
up to ~ 18% of cases in 2015, and up to ~ 25% in 2016. Ob¬ 
serving for 4 days (600 fields) this becomes 25% in 2015 and 
~ 32% in 2016. These numbers assume that nearby sGRBs 
will be comparible in luminosity to those detected to date, 
which is unlikely therefore they represent upper limits on 
the fraction we can identify. 

In conclusion: the X-ray wavelength offers a promising 
means to identify the EM counterpart to a GW trigger, due 
to the low rate of bright transients compared to other wave¬ 
lengths, and the fact that X-ray facilities are in orbit, and 
thus not subject to problems of daytime. Swift, with its flex¬ 
ible planning and rapid slewing abilities, is able to cover sig¬ 
nificant portions of the sky, and by conducting a series of 
short (50-500 s) exposures, has a reasonable chance of iden¬ 
tifying any X-ray counterpart to the GW trigger that may 
exist. 


APPENDIX A: THE XRT SIMULATIONS 

To simulate the actual follow-up observations, we used the 
following approach (when we describe a field as ‘observable’ 
below, we mean that it is observable by Swift-XRT for at 
least e s): 

(i) Set a counter 1 = 1. This keeps track of how many 
iterations over all N fields have been attempted. Set the 
time to Td (Section 121). 

(ii) Identify the most probable field observable. If no field 
is observable, identify which one will become so soonest, and 
wait until then. 

(iii) Increase the exposure of the current field by e s. 

(iv) For each unrelated source in the field which has not 
yet been detected, determine the probability that the source 
is now detected, based on the its count-rate, the exposure in 
this field so far, and the completeness of the detection system 
llEvans et al. I I 2 OI 4 section 7 and fig. 14). Then generate a 
random number between 0 and 1 to determine whether the 
source is detected. 

(v) If the GRB lies in the field, determine how many 
counts are detected from it in this exposure, from the 
redshift-corrected afterglow template defined in Section [S] 
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Add this to any counts detected in previous (or overlap¬ 
ping) observations of the GRB. Then determine, whether or 
not it is detected, as for the unrelated sources. If it is de¬ 
tected, determine the mean count-rate of the detection, and 
compare this to the RASS upper limit for this field. Also, 
construct a light curve of the GRB, and determine whether 
fading can be determined with 3-a confidence. If the GRB 
is above the RASS limit, or fading, mark it as ‘identified’ 
and record the time. 

(vi) If all fields that are not too close to the Sun, Moon or 
orbital pole to observe, have been observed / times, increase 
I. 


(vii) If there are any fields which overlap the current field, 
are observable, and have been observed < I times, identify 
the one with the greatest probability, and select it. 

(viii) If there are no overlapping fields matching the above 
criteria, identify all currently observable fields that have 
been observed fewer than I times, and select the one with 
the highest probability. If any field which has been observed 
fewer than 7 — 1 is observable (e.g. it was Sun constrained 
during the earlier rounds of observations, but is now ob¬ 
servable), select this field in preference to the other fields, 
regardless of probability. 

(ix) If no fields are currently observable, determine when 
the next field that has been observed fewer than I times 
will become visible, and wait until this point, then select 
that field. 

(x) If the time is less than 15 days since the trigger, go 


to step (iii) otherwise, end the simulation. 


Step (vii) exists for efficiency reasons: by observing over¬ 


lapping fields consecutively by preference, the ratio of time 
spent observing to time spent slewing is maximised. The 
GW probability maps have smooth contours, therefore the 
only time this step results in selecting a field with lower 
probability than a non-overlapping field is when the map 


bifurcates, in which case not including step (vii) would tend 


to result in Swift constantly slewing between the two dis¬ 
jointed regions, wasting much valuable time. 

The simulations assume that to move from a field to an 
overlapping one requires 25 s of ‘dead time’ between the end 
of one observation and the beginning of the next. This value 
was determined by examining data from the S wift tiled ob¬ 
serva tions used to respond to neutrino triggers jEvans et al.l 
[ 201 ^. The situation when performing a slew between dis¬ 
jointed fields is more complex, as it depends on various de¬ 
tails of the slew, which is beyond the scope of this work 
to simulate. Based on an analysis of historical Swift obser¬ 
vations, we assumed a slew rate of 0.77° s“^, plus 30 s of 
overhead. 

After each observation, the simulation records the time 
since the trigger, how many unrelated sources have been de¬ 
tected so far, and how many of those are not in the RASS, 
and whether or not the GRB has been detected, and iden¬ 
tified. 
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